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A binary partition of a positive integer n is a partition of n in which each part 
has size a power of two. Let b(n) denote the number of binary partitions of n. Since 
a binary partition either has a part of size 1 or else is twice a binary partition of 
n/2, b(n) satisfies the recurrence 

bin) = bin — 1), n odd; 

(1) 

b(n) = b(n — 1) + b(n/2), n even. 

The generating function J] b{n)x n = 1/ 11(1 — x 2 ) and small values of b(n) were 
written down by Euler [2j §50], and Mahler [3j and de Bruijn pQ gave increasingly 
good asymptotics for b(n) and for partitions into powers of r other than 2; this is 
sometimes called "Mahler's partition problem." 

In this note we first use a variation on the recurrence |T]) to construct a Gray 
sequence on the set of binary partitions themselves. This is an ordering of the set 
of binary partitions of each n (or of all n) such that adjacent partitions differ by 
one of a small set of elementary transformations; here the allowed transformatios 
are replacing 2 k + 2 k by 2 fc+1 or vice versa (or addition of a new +1). Next we give 
a purely local condition for finding the successor of any partition in this sequence; 
the rule is so simple that successive transitions can be performed in constant time. 
Finally we show how to compute directly the bijection between k and the fcth term 
in the sequence. 

Thanks to Donald Knuth for requesting Theorem Q] and to Richard Stanley for 
responding to an early draft by asking about the possibility of Theorem [3l 

1. Construction of Gray sequences 

Theorem 1. For each positive integer n, the binary partitions of n can be arranged 
in a sequence B(n) such that adjacent partitions differ by an operation of the form 

■■■ + 2 k + 2 k + ■■■ < — > ■■■ + 2 k+1 + ■■■ . 

Moreover the sequence runs first through all partitions ending with ■■■ + !, and next 
through all partitions ending with ■ ■ ■ + 2. 

Proof. We give a recursive construction for B(n). We let Q(n) and S(n) denote 
the first and last partitions in the sequence (reserving R for later use); Q{n) will 
always be 1 + 1 + • • • + 1. 
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Here x 2 indicates doubling each part of a partition. Note that in the 2 mod 4 case, 
the f sequence appears in reverse order. 

The same logic that justified formula (fT]) for b(n) shows that the above sequences 
certainly contain all binary partitions. By induction, we only need to check that 
the two points of concatenation obey the adjacency condition. 

For n = mod 4, S(n - 1) + 1 = S(n - 2) + 1 + 1 = Q(f - 1) x 2 + 1 + 1 (since 
n — 2 = 2 mod 4), which is 2 + • • • + 2 + 1 + 1, since Q(m) is the all l's partition for 
all m. Meanwhile, Q(f ) x2is2 + -- - + 2 + 2, and is thus connected to S(n — 1) + 1 
by a 1 + 1 2 move. 

For n = 2 mod 4, S(n - 1) + 1 is again 5*(n - 2) + 1 + 1, but now n - 2 = mod 4, 
so this is S(% - 1) x 2 + 1 + 1. On the other hand, 5(f) x 2 is S(f - 1) x 2 + 2 
(because f is odd), so S(n — 1) + 1 and S(§) x 2 are also connected by a 1 + 1 «-> 2 
move, and we are done. □ 

As n increases, the head of B{n) remains unchanged aside from adding +ls to 
each partition. So there is a single infinite sequence B = B\, B2, ■ ■ •, beginning 

0, 2, 22, 4, 42, 222, 2222, 422, 44, 8, 82, 442, 4222, 22222, 222222, 42222, 4422, 822, 84, 444, . . . 

such that each B(n) is just the initial substring of B of partitions summing to < n, 
padded with the appropriate number of Is. 

B is a list of all binary partitions with only even parts, so halving each one gives 
a Gray sequence 6/2 of the binary partitions of all n. Here the notion of a legal 
transition must be expanded to include the operation P — > P + 1, the remnant of 
the transition 1 + 1 — ^ 2 after dropping all Is and then halving. By the construction, 
the subsequence of partitions in B/2 with constant sum n is identical to B(n) above 
if n is even, and is the reverse of B(n) if n is odd. 

2. Stepping through the sequence 

Given a partition P in the sequence B, we can calculate the partition which 
comes before or after P easily. We will give explicit maps 0+ and 0_ which take 
a binary partition and return which of the rules 2 k + 2 k ^ 2 k+1 transforms P into 
its successor or predecessor in B. 

Looking back at the construction of B from the B(n), we see that for most P, 
the adjacent partitions <fi±(P) are the same size as P; by "size" or \P\ we mean 
the sum of the parts of P when viewed as an element of B, so after discarding any 
parts of size 1. The only exception is when we try to apply to the last partition 
of size n, which we called S(n), or to apply 0_ to the first partition in B of size 
n, which we now name R(n). In the construction of B{n) this is the first partition 
which results from a x 2 operation. 

From the construction of the B(n), we can calculate: 

2 + ■•■ + 2, 
2 + ••■ + 2, 

S{n/2) x2 = 5(n/4) x 4 = ... 
2° H + 2° , where n = 2 a b with b odd, 

b 

S{n-2) + 2. 

Conveniently, S(n) for n = 2 mod 4 fits the mod 4 pattern as well. 



(3) 



For n = mod 4, R(n) 
For n = 2 mod 4, S(n) 
For n = mod 4, S(n) 



For n = 2 mod 4, R(n) = 
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So the special cases are easily identified, and we already know which transfor- 
mation rules to apply to them: 

MS(n)) = 1 + 1-2 
y > <f>-{R(n)) = 2—1 + 1 

Of course, all parts of size 1 are suppressed in B, so the effect of the rule 1 + 1 — 2 
is that a part of size 2 appears from nowhere, and the apparent size of the partition 
increases. 

Aside from these special cases, all other transformations can be determined re- 
cursively: 

f 2x<MLP/2j), |P|=0mod4 
W <P±^)- \ 2x<^(LP/2j), |P| = 2mod4 

Here |_ - J denotes deleting all parts of size 1 to again get an even binary partition, 
and the =F in the 2 mod 4 case accounts for the reversal of S(f ) in the definition 
of B(n). The multiplication by 2 acts on the rule returned, so 2 x (1 + 1 — 2) is 
2 + 2 — 4. 

Unravelling the recursion leads to a startlingly quick way to step through B. 

Theorem 2. Suppose you are given a binary partition P written . . . d 3 d 2 did , 
where the the digit dk is the number of parts of P of size 2 fc . Since we are working 
in B, we will ignore the value of the Is place d , except to assume it is at least two 
if the transformation 1 + 1 — 2 is needed. 

• Let i be maximal with di > (i = if P = %). 

• Let j be the second largest integer with dj > (or j = if none). 

• Let e be ( — 1)^-- dfe summing over 1 < k < i — 1. 

Then the transformations <j>± act on P according to the following rules: 
If (di, dj, e) is. . . then <p± says. . . 

(a) ( 1 , *,+l) dil,di-\X\, split a largest part 

(b) any other (odd,*, =pl) di + i |, rf, H, merge two largest parts 

(c) (odd, 1,±1) dj i, dj-i ft; split a 2 nd -largest part 

(d) any other (odd,*, ±1) dj + \ |, dj JJ,, merge two 2 nd -largest parts 

(e) (even, *, +1) di j, di—i ft, split a largest part 

(f) (even, *,±1) di+it,rfjj|, merge two largest parts 

Here |, I indicate an increment or decrement by one, and ft, II &?/ ^wo. 

Note in particular that every transition in £> involves merging or splitting one of 
the largest two sizes of parts in P\ As a result, successive updates can be done in 
place in constant time, in a well-chosen data type (where you keep track of e, and 
you never need to search for dj, e.g. a linked list of nonzero dk-) 

Proof. The recursive part of the definition of (f>± is handled trivially in this notation. 
The operation [P\ just shifts each digit to the right and forgets do, and the value 
of \P\ mod 4 depends only on the parity of the 2s digit d\. So we keep dropping 
rightmost digits until we arrive at one of the base cases 4> + (S(n)) or 
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and just remember e to know if we've reversed direction an even or odd number of 
timefl 

The rest of the proof consists of identifying which of the base cases listed in 
equations ([3]) and (fJJ is the destination of each partition: 

• Rule (d) is for partitions that end at (j> + (S(n)). The (odd) lead digit di is 
the (odd) number of identical parts in S(n) = 2 a + ■ ■ ■ + 2 a . We reach the 
base case as soon as the second-largest nonzero digit dj is deleted, and the 
transformation 1 + 1 — > 2 therefore joins two 2 J s into a 2 J+1 . 

• Rule (c) takes care of the exception where we run into <f>-(R(2n + 2)) one 
step before we would otherwise reach cj) + (S(n)). 

• Rule (e) is for partitions that end at <fi-(R(n)) when n = mod 4, an 
even- length sum 2 + • • • + 2. 

• Rule (a) is for the somewhat special case <fi-(R(2)). R(2) = 2, the only 
time R{n) is a partition with an odd number of parts, since S(0) = 0. 

• Rules (b) and (f) are the fall-through cases: the recursion avoids all base 
cases until it gets to </>+(0), to which the transformation 1 + 1 — > 2 is applied. 

□ 

For example, let us compute the next several partitions in B beginning with 



256 5 32 2 16 1 4 


4 2 3 


(so 


1^1 


= 1382). 
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(d) 
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(b) 
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(a) 
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Note that e changes sign under rules (b) and (c), and under (d) unless i = j + 1. 

3. Calculating individual terms 

The recursive definition leading to the sequence B = B\,Bi,... allows us to 
explicitly compute the bijection k «-> Bk- First we introduce an alternate notation 
for even binary partitions. 

Definition. Let P be an even binary partition (that is, with no parts of size 1). 
Define the trail of P, t(P) = r = tq, t%, T2, ■ ■ ■ by 

n = \[P/r\\ 

where again |_-J indicates deleting any parts of size 1 or smaller to again get an 
even binary partition. 

In other words, is the size of the partition after i iterations of the map "halve 
all parts and delete parts of size 1." For example, if P = 88422, then [P/2\ = 442 
and L-P/4J = 22, so t(P) = 24, 10, 4, 0, . . ., and we henceforth omit the trailing 0s. 



The fact that we ignore do in the definition of e perhaps reflects a moral imperfection in our 
definition of B{n): maybe for odd n it should be the reverse of B(n — 1), which then eliminates 
the reversal in the 2 mod 4 case. 
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The partition P is easily recovered from t(P): the number of parts of P of size 2 1 
is Tj-i/2 — Ti, the number of parts of size 1 dropped on the ith iteration. 

Theorem 3. Given just the integer k, the trail of Bj~ can be determined as follows: 

(i) \Bk\ is the smallest n such that k < b(n). 

(ii) With n = \Bk\ as above, [Bk/2\ = Be, where 

e = k-b(n-2), n = 0mod4, 
I = b(n) + l-k, n = 2 mod 4. 

Conversely, the trail r = tq, t\,T2, . . . corresponds to the kth partition Bk if its 
truncation 7i, T2, T3, . . . corresponds to the ith, where 

k = 6(t -2)+£ r = mod 4, 
k=b(T ) + l-t, t = 2 mod 4. 

This recurrence and the base case t = 0, 0, . . . corresponding to B^) = B\ = 
suffices to determine the location in B of the binary partition with any given trail. 

For example, partition 88422 with trail 24,10,4 appears at position 6(22) + 
6(10) + 1 - (6(2) + 6(0)) = 86. And to find 23i 23 456789, we calculate that 

123456789 G (6(646), 6(648)] note 648 = mod 4, 
123456789 - 6(646) G (6(304), 6(306)] note 306 = 2 mod 4, 
6(306) + 1 - (123456789- 6(646)) G (6(120), 6(122)] etc, 

and its trail is 648, 306, 122, 58, 28, 14, so the partition is 64 7 32°16 1 8 3 4 31 2 18 . This 
amounts to a writing of 123456789 in terms of values of b(n) with n = 2 mod 4: 

6(646) + 6(306) + 1 - (6(122) + 1 - (&(58) + 1 - (6(26) + 6(14) + 1 - (&( ))))) 

It is possible to work with this representation directly, using sign changes to track 
the mod 4 behavior, but statements tend to be inelegant. 

Both maps rely heavily on the values of the function b(n) which counts the 
number of binary partitions of n, and for which we have no closed form. From 
a computational complexity point of view this is inevitable, since from the map 
k — > Bk the values of b(n) can be determined easily. In practice the computation 
of b(n) by recurrence |T]) is inexpensive. 

Proof. Since the trail of any even binary partition P begins with tq = \P\, the given 
maps k — > r and r — > fc are clearly inverse, so it suffices to show either direction. 

The map fc — > r(Bk) follows directly from the construction of the sequences B(n). 
Suppose n is even. After deleting all parts of size 1 from the partitions in all B(n), 
the first b(n — 2) terms of B(n) are exactly B(n — 2), and the remaining 6(n/2) 
terms are the even binary partitions of n. So these partitions occur in B as a solid 
block of terms with indices in (b(n — 2), 6(n)], justifying (i). 

For (ii), recall that each of the terms in this block was obtained by doubling 
some term in B(n/2); the [-J operation deletes all parts of size 1 and allows us to 
work in B instead. If n = mod 4 then we doubled Bi to get the £th term in the 
block, £>b( Tl _2)+^ 5 while if n = 2 mod 4 then B(n/2) was reversed, and we doubled 
Be to get the ith term from the end, Bu n )+x-l- ^ 
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